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The dynamics of a tracer molecule near a fluid membrane is investigated, with particular emphasis 
given to the interplay between the instantaneous position of the particle and membrane fluctuations. 
It is found that hydrodynamic interactions creates memory effects in the diffusion process. The 
random motion of the particle is then shown to cross over from a "bulk" to a "surface" diffusive 
mode, in a way that crucially depends on the elastic properties of the interface. 
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I. INTRODUCTION 



A largely unsolved problem in soft materials is to un- 
derstand how the random motion of a colloid is corre- 
lated with the viscoelastic properties of the embedding 
fluid 0. Since complex fluids store elastic energy, the 
instantaneous position and velocity of the Brownian par- 
ticle depend on its prior history. As a signature of this 
"memory" effect, the friction force experienced by the 
particle is, in ageneralized Langevin description, non- 
local in time @>0|. To illustrate this point, we consider a 
system simple enough that allows us to extract the mem- 
ory kernel. We focus in this paper on the diffusion of 
nano-particles near a liquid-like membrane. Fluid mem- 
branes are soft surfaces, self-assembled from surfactant 
solutions. In living systems, they divide the cell into din- 
stinct compartments and incorporate a large amount of 
proteins involved in signalling activities. Interactions be- 
tween nano-sized particles and fluid membranes are ubiq- 
uitous in phenomena such as antibiotic delivery, nano- 
particle toxicity or virus entry in cells |^. This issue is 
relevant not only from biological or biophysical stand- 
points, but also in the formulation of many cosmetics 
and pharmaceutical products 5]. At a more fundamen- 
tal level, recent experiments on the dynamics of particles 
near a membrane have revealed peculiar kinetics effects 
such as anomalous diffusion [(| , or hierarchical transport 
mechanisms Besides, anisotropic coherent motion of 
particles induced by an alterning electric field has also 
been shown to enhance the smectic order in a lamellar 
phase of membrane Q . From all these perspectives, a de- 
tailed investigation of the Brownian motion in the vicin- 
ity of a membrane appears to be essential. 

Interfacial contributions to transport coefficients have 
been evaluated long ago for elementary geometries 0. 
Far from the surface, the mobility of a spherical particle 
of radius a is /io = 1/ (6^770) , with 77 the shear viscosity 
of the embedding (simple) fluid. When a horizontal wall 
is located at distance zq from the particle, the fluid flow 



resulting from the motion of the bead "reflects" from 
the interface back to the particle, thus exerting an ad- 
ditional friction force. At leading order, hydrodynamic 
corrections lead to a position-dependent mobility tensor 
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and fj,ij = for i =/= j. Note that within the assumption 
of point-like particle a <C zp , th e exact shape of the col- 
loid is essentially irrelevant ^(j- I n Eqs. ljTa|) and Qlbfl. 
particle mobility is affected in a manner that depends 
on the nature of the interface. For a solid surface, the 
no-slip boundary condition leads to |j| 
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During the last decade, precise experimental measure- 
ments of particle diffusivity near solid surfaces have 
shown excellent agreement with theory 0, 0, 0, 0] . 
If, however, the interface separates two liquids having the 
same viscosity 77, the corrections are 
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More recently, the effects of partial slip [lfj] as well as 
fluid compressibility |17| have also been included in the 
calculation of the mobility tensor. 

Unlike surfaces having a frozen structure, fluid mem- 
branes are responsive materials. The Brownian motion 
of a particle near a membrane is then expected to be dy- 
namically coupled to the conformations of the interface. 
In Section ^ we consider the problem of free diffusion 
near a membrane, with particular emphasis given to hy- 
drodynamic interactions. We show in Section II I II that 
the instantaneous solution of the creeping flow equations 
depends on the whole history of both colloid's and in- 
terface's motions. Consequences for the diffusion process 
are then discussed in Section llVI In particular, we find 
that the random motion of a free colloid crosses over from 
a "bulk" to a "surface" diffusive mode. Finally, we come 
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back to the relationship with experiments and draw some 
concluding remarks in Section|Vl the technical details be- 
ing relegated to the appendices. 



II. HYDRODYNAMICS NEAR A FLUID 
MEMBRANE 

A. Physics of membranes 

Fluid membranes belong to the general class of soft 
objects (including polymers, gels . . . ), that deform easily 
when submitted to external stresses, and undergo ther- 
mally excited shape fluctuations to increase their config- 
urational entropy. At length-scales much larger than the 
bilayer thickness (i.e., micrometers vs. nanometers), a 
fluid membrane can be described as an infinitely thin liq- 
uid layer with essentially no in-plane shear viscosity |18| . 
The instantaneous conformation of the almost flat mem- 
brane is specified through the displacement field h(p,t), 
with p = (x.y). For our purpose, it is more convenient 
to use the (2D) Fourier representation 



h q (t) = J d 2 pexp[-iq.p]h(p,t) 
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with q = (q Xl q y ). As first explained by Helffich, the large 
variety of shapes and topologies assume by membranes is 
governed by bending energy. For relatively smooth defor- 
mations, the statistical mechanics of a single membrane 
is based on the following Hamiltonian [l^ 
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with k the bending rigidity. In this description, the mem- 
brane is confined in a harmonic potential. The in-plane 
correlation length £|| plays the role of an external control 
parameter setting the mean square displacement of the 
membrane: (h(p) 2 ) 1 / 2 cx [k B T / n) 1 / 2 ^, with T the ab- 
solute temperature and k B the Boltzmann constant. As 
recalled in Appendix IaI the relaxation rate of a fluctu- 
ation with wavevector q follows the dispersion relation 
Lu q = K(q 4 + £|^ 4 )/(477g) 20]. Overdamped surface waves 
are further characterized by their longest relaxation time 



4 v 



(6) 



For a typical bilayer in water with 77 — 10 3 Pa.s and 
£|| = 0.5 /xm, this relaxation time varies between r m ~ 
10~ 3 s for k — 4.10" 19 J and r m w 0.1 s for k = 4.10" 21 J. 

In order to describe the dynamic organization of a col- 
loidal suspension near a membrane, r m has to be com- 
pared with the typical time needed for a particle to 
diffuse over its own radius 
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This time-scale ranges from r<2 w 10~ 6 s for a particle 
with radius of a few nanometers, up to t<j w 1 s for mi- 
crometer beads. Comparison of both time-scales defines 
two diffusion regimes in a straightforward manner. When 
Td/T m 3> 1, i.e. when a 3> £|| or k 3> fcsT, the membrane 
appears essentially flat to the colloid and the system is 
dominated by the relaxation dynamics of the elastic in- 
terface. We call this regime the weak fluctuation (WF) 
regime. On the other hand when Td/r m <C 1, the bead is 
strongly advected by the random flow caused by thermal 
undulations of the membrane. We will refer to the strong 
fluctuation (SF) regime, that will be the particular scope 
of investigation of this paper. 



B. Linear hydrodynamic 

For small-amplitude motions, the flow velocity u and 
the local pressure p are governed by the linearized Navier- 
Stokes equation 



?yV 2 u -Vp + f + * = , 
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together with the incompressibility condition V.u = 0. 
Here, f is an external force density causing the fluid mo- 
tion, and <& the restoring force density exerted by the 
deformed interface. The latter is always perpendicular 
to the membrane. Within the assumption of smooth de- 
formations, the force density is given in Fourier space 
by the functional derivative <fr = (0, 0, —SH m /Sh*). The 
creeping flow Eq. © is then solved assuming that the 
membrane is not permeable to the fluid, a condition quite 
sensible on experimental time-scales (minutes or hours). 
Explicitely, this requirement reads 



dh 
~dt 
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The hydrodynamic problem defined in this way is both 
time-dependent and highly non-linear due to the fact that 
h(p,t) is unknown. Althought it cannot be solved ex- 
actly (excepted by numerical methods), approximate so- 
lutions can be found. For instance, iterative methods 
have proved to be fruitful when the deformation of the 
interface is asymptotically small 21). The zeroth-order 
term usually corresponds to the condition u z (p, 0, t) = 0, 
thus representing the motion of a spherical bead near a 
flat surface. The ensuing velocities and stresses can then 
be used to determine a first nonzero approximation for 
the deformation of the interface [2]]]. This strategy is 
however limited as it only accounts for the first "image" 
correction to hydrodynamic interactions. Instead, we fol- 
low a different route and assume that the fluid motion 
caused by thermal fluctuations of the membrane is only 
slightly perturbed by the presence of the colloid. The 
dynamics of the surface enters the problem through the 
condition 
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The consistency of this approach is then ensured as long 
as non-linear contributions in the deformation field can 
be neglected. In fact, this condition has the same range 
of validity as the harmonic description for the membrane 
Hamiltonian l22l. 



colloid, and so forth. The infinite sum (u> — iujq)^ 1 — 
i/u) q X)^Lo(~*)" ( w / ' w <?)" ^ s * ne expression of this infinite 
series of "reflexions" of the original point force on the 
interface. 



III. MOBILITY OF A SPHERE NEAR A FLUID 
MEMBRANE 

A. Green's function 

Standard Green's function methods are then used to 
solve Eq. JHJ. We first evaluate the response function, 
G, to a point-like force fo<5 (r — r'). The fluid flow corre- 
sponding to a given force field f(r) is then obtained from 
the convolution 



u(r,t) = Jdtj d 3 r'g(r,r',t-t'){ (r',i') 
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For the hydrodynamic problem defined above, the Os- 
een tensor splits into two contributions, Q = + AG- 
The first term, G^\ is the well-known free space Green's 
function and is given in Appendix ^ The second contri- 
bution, AQ, describes the fluid flow caused by the elastic 
response of the membrane. While the problem is trans- 
lationaly invariant in time and in the directions parallel 
to the surface, AG is expected to depend on the perpen- 
dicular coordinates z and z' separately. Incidentally, we 
remark that the boundary condition Eq. (|I0|l is more con- 
viniently enforced if we use the set of variables {q, z, uj}, 
where we define the (time) Fourier transform of any func- 
tion F(t) as F(uj) = dte~ lut F{t). The main lines 
of the derivation are given in Appendix iBl We find after 
some algebra 



AGm (q, z, = 



-. — lk(q,z)"/i{q,z')Mki , 

Arjq ui — iuj q 

where the matrix A4 has diagonal elements Aikk — 1, and 
off-diagonal terms M xy = M yx = 1, M xz = M yz = —i, 
and M.zx — M-zy = i. The functions 7^ are given by 
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Remark that these expressions are consistent with 
the symmetry property of the Green's functions: in- 
deed, it can be shown from general arguments that 
Gki (r, •*') = Gik (r', r) j2^|. This is re-expressed in terms 
of our particular choice of variables as Gki (q, z, z' , to) — 
Gik {— q, z 1 ', z, lo), property that can be checked easily. 

In Eq. l(T21) . the prefactor (iv — iujq)^ 1 is a clear sig- 
nature of hydrodynamic scattering effects. Indeed, the 
fluid flow resulting from a displacement of the parti- 
cle exerts stresses that deform the membrane. Relax- 
ing back to its equilibrium position, the membrane cre- 
ates a backflow that in turn perturbs the motion of the 



B. Mobility tensor 

The next step in the derivation consists in identifying 
the mobility tensor of a sphere. To this aim, we still have 
to enforce the no-slip boundary condition for the fluid 
flow on the surface of the colloid. If we note v(ro) and 
f2 respectively the translational and rotational velocity of 
the Brownian particle, ro being the position of its center- 
of-mass, then the fluid velocity is 



u(r + a) = v(r ) + J7 x a 



(14) 



for any vector a scanning the bead's surface. The friction 
force Fh exerted by the liquid on the particle is then 
obtained from 



F H = - J dVi (r) = - I dSt(a) 



(15) 



where the integral is running over the surface of the 
bead. In this expression, f(r) = f(a)5(|r — ro| — a) is 
the force density exerted by the particle on the fluid. 
Integrating Eq. (fT4*jl over the particle's surface, the an- 
gular velocity cancels out and one obtains, together 
with Eq. (|11|) . a linear relation between the friction 
force and the velocity of the particle. This relation de- 
fines the (frequency-dependent) mobility tensor, v(u>) = 
— J1(lu)Fh (tt>). The mobility can be written as the sum 
of two terms, #jm(zo,w) = ^ 4; + AJl k i(z ,uj), with 
/iQ = l/(67r?ya) the bulk value. In the limit of a point- like 
particle a zq, we obtain at leading order 



d 2 q 
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From Eqs. I|13a(l - l|13c|l . one can easily convince oneself 
that there are no cross-contributions: AJl^i = if k ^ I. 



1. Static mobility 

The integral Eq. ||TSJ| can be performed easily when 
u> = 0. As can be noticed from Eq. I|12|l . the static mo- 
bility does not depend on the elastic properties of the 
interface. Indeed, we find 



Hxx(zo,0) = fiyy(zo,0) = Mo U 
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We therefore recover the mobility of a sphere near a flat, 
liquid interface given in Eq. Q- 
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2. Low-frequency limit 

For finite frequency, AJlki is a complicated function of 
the reduced variables wr m , a/zo and £\\/zq. At this point 
of the discussion, it is usefull to introduce the friction 
coefficient defined as the inverse of the mobility tensor, 
Cki = (/^ -1 )fcz- Obviously, the friction tensor is also di- 
agonal. In the low- frequency limit u)T m <C 1, the matrix 
elements can be expanded and we find, at lowest order, 



Ckl (zo,w) = Co 



1 



(18) 



with Co = I/Mo- The prefactors flu depend on the dis- 
tance to the membrane. We get f3» = and f3± — 
37ra/(16^||) when £y 3> z§, whereas in the other limit 
£|| <C zo we obtain (3u — 9o^||/(64zq) and (3j_ = 
27<||/(32z 2 ). 

Writing the complex friction as C( w ) = 67ny(w)a, it 
appears that Eq. I|18|l describes a Maxwell fluid with 
anisotropic, frequency-dependent shear viscosity Q]. The 
corresponding relaxation time is essentially set by the 
membrane relaxation time r m , multiplied by a geomet- 
ric prefactor that depends on the distance to the elastic 
surface. This relaxation time arises from the delay in the 
response of the membrane to a deformation caused by 
the fluid flow. 



IV. DIFFUSION NEAR A FLUID MEMBRANE 

We are now in a position to study the diffusion of a 
tracer particle in the vicinity of a membrane. Because of 
hydrodynamic interactions, the friction force experienced 
by the Brownian particle is non-local in time. The cor- 
responding generalized Langevin equation for a particle 
of mass to reads 



dv 
i— 

dt 



dt'({t - t')v(t') + F(t) , (19) 



with F(i) the noise term. All the statistical properties of 
this random walk can, in principle, be calculated given 
the probability distribution of the random force. How- 
ever, the problem is still involved because of the non- 
linear z-dependence of the friction term. To avoid this 
complication, we assume that the particle remains close 
to its original position during the time of observation — 
the validity of this assumption being discussed later. Ac- 
cording to the fluctuation-dissipation theorem, the ran- 
dom force F(t) can be choosen with zero mean value and 
correlations given by [24l |25| 



F k (u)Fi{u')) = 2k B TRe[Jl kl (uj)} x 2ttS(lj+u') . (20) 



We first focus on the motion in the direction perpen- 
dicular to the surface. In the overdamped limit, i.e. for 
time-scales much longer than r = C,o/m, inertial terms 




Reduced time 

FIG. 1: Excess mean square displacement A(t)/A(t — > oo) 
(see text) as a function of the reduced time t/r m - The curves 
correspond to different values of the reduced distance d = 



can be neglected. Defining the mean square displacement 
(MSD) as (Sz 2 ) = ((z(t) - z ) 2 ), with z = z(t = 0), we 
obtain 



(Sz 2 )=2D a t- kBT 1 ' 2 " 1 
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dqe- 2qZ0 (1 + qz f 



X it + 



(21) 



Here, Do = l^oksT is the bulk diffusion coefficient. This 
result presents some interesting features. At short times, 
the MSD behaves like 



{Sz 2 ) ~ 2D t , 



(22) 



for t <C t to . In this limit, the particle does not "feel" 
the membrane. On the other hand, one retrieves for 
long time-scales the mobility of a particle near a non- 
deformable liquid interface 



(5z 2 ) ~ 2D ( 1 
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(23) 



for t 3> r m . In this limit, the MSD is thus independent 
of the elastic properties of the surface 1 . Both the unex- 
pected short-time and the long-time behaviours can how- 
ever be explained from the definition of r TO . As can be 
noticed from Eq. © , the limit t ^> r m actually coincides 
with the limit of an infinitely rigid interface K — » oo. 
Therefore, one ultimately expects to recover the result 
for a flat interface when t — ► oo. On the other hand, the 
limit t <C T m corresponds to a membrane with vanishing 
bending rigidity k — * 0. One can easily check that all 
the corrections discussed so far actually die out in that 



1 In fact, the long-time diffusion coefficient could also have been 
obtained from the general relation D± = kBTJl±(Lj = 0), com- 
pare with Eq. 117bl . 
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limit, so that we recover the bulk diffusion coefficient 
when t — ► 0. 

For intermediate time-scales, one crosses over continu- 
ously from the fast to the slow diffusion mode. Defining 
A(t) = {8z 2 ) — 2D±t, we compute numerically and plot in 
Fig. [2 the excess MSD A(t)/A(t -» oo) for different val- 
ues of the ratio d= zq/£h. When the particle is far away 
from the membrane, surface deformations and fluctua- 
tions are essentially irrelevant and the particle reaches 
the "long-time" diffusive regime on a time-scale much 
shorter than r m . However, as the particle approaches 
the membrane, the strength of the hydro dynamic cou- 
pling between the motion of the particle and the fluctu- 
ations of the membrane increases. The effect of the fric- 
tion kernel on the Brownian motion of the colloids can be 
seen on much longer time-scales, and the cross-over time 
increases up to ~ r rn for zo ~ £|| . As expected, hydro- 
dynamic effects are all the more relevant as the bending 
rigidity (or the confinement potential) of the membrane 
are small, corresponding to large values of r m . 

As far as the x and y directions are concerned, one just 
has to replace the factor (1 + qzo) 2 by qz /2 in Eq. i(2T|) 
to get the corresponding MSD. The results are essentially 
the same as for the z direction. In particular, one recovers 
the diffusion coefficient Du = fcsT^ (l — 3a/(32z )) in 
the long-time limit for the same reasons as stated above. 

Finally, we still have to comment on our approximation 
that neglects the z-dependence of the friction coefficient. 
This assumption is valid as long as (6z) <C Zq, that is for 
t <C r = azQ7]/(kBT). Therefore, the results described 
above can be observed provided that To > T m , so that we 
reach the condition 



ksT az, 



(24) 



For an initial height zq = 10£ii and a typical bending 
rigidity k = lOfc^T, this condition is readily satisfied. 



DISCUSSION AND CONCLUSION 



wall, whereas the long time-scale correlations come from 
fluctuations of the membrane itself. The latter create in 
turn intensity fluctuations by modulating the number of 
detected particles. The results are interpreted by consid- 
ering that diffusion is anomalous due to the collisions of 
the particles with the membrane l2fij . This mechnism 
is however not incompatible with our results: indeed, col- 
lisions have been neglected so far as they are not easily 
included within a Langevin formalism. In order to get a 
more accurate description of the system, both collisions 
and hydrodynamic interactions should therefore be taken 
into account. Note also that the hydrodynamic coupling 
described in this work is expected to be most pronounced 
for a ~ zq ~ £|| . A more carefull study, presumably nu- 
merical, of these various points is therefore required. 

In conclusion, we have investigated the motion of a 
Brownian particle in the vicinity of a liquid-like mem- 
brane. The formalism developped in this work allowed 
us to accurately account for membrane- induced hydrody- 
namic interactions. Obviously, this study might be rele- 
vant from a biological perspective as diffusive properties 
of nano-particles near a membrane have been shown to 
crucially depend on the elastic properties — or equiva- 
lently the composition — of the bilayer. In particular, it 
would be interesting to see whether the predicted mem- 
ory effects can influence the various kinetic processes tak- 
ing place near a membrane. 



APPENDIX A: FLUCTUATING 
HYDRODYNAMICS NEAR A FLUID 
MEMBRANE 

Thermal undulations of the membrane actually origi- 
nate from hydrodynamic fluctuations of the embedding 
fluid. In this appendix, we re-derive the well-known 
equation of motion for the membrane [Tfll l2f)j | directly 
from fluctuating hydrodynamic. In Fourier representa- 
tion, the creeping flow equation satisfied by the flow ve- 
locity u(k,i) = / d 3 re- ik - r u(r,t) is [HHg 



To summarize, we have shown that the random motion 
of a colloid near a soft membrane is a non-Markovian 
process. Indeed, the delay in the response of the elastic 
interface to hydrodynamic stresses induce memory effects 
that are relevant over almost three decades in time (see 
Fig.nj. Consequently, the Brownian motion of a colloidal 
particle might locally appear anomalous with MSD given 
by (5r 2 ) ~ V . But as explained above, hydrodynamic in- 
teractions are not expected to lead to "true" subdiffusive 
or superdiffusive behaviours since the MSD is always lin- 
ear in time when t — > oo. 

Actually, anomalous diffusion has been observed in re- 
cent fluorescence correlation spectroscopy experiments 
close to a vesicle 0. It has been shown that the in- 
tensity autocorrelation function is the sum of two contri- 
butions. The short time-scale correlations arise from the 
(regular) diffusion of the particles near an impenetrable 



u(M) 
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x-k£)#(k,t) + g(M), (Ai) 



with X the identity tensor, kk|y- = k%kj the dyadic 
product, and k = k/A;. The restoring force density 
«1? = (0, 0, <& z ) is given, in Fourier representation, by 
<&z(k, u) = —SH m /Sh* = —E q h q (uj), with the energy 
density E q = n(q 4 + £,7) and q = (k x ,k y ). In Eq. (23) j 
the field g(k, t) is a stochastic variable characterizing 
thermal fluctuations in the liquid. Following Ref. [28| . 
it can be shown that the random noise has a Gaussian 
distribution with zero mean value, (g(k, t)) = 0, and 
correlations 



< 5i (k,%(k',o> 
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where the square brackets denote average over an equi- 
librium ensemble. 

In order to derive the equation of motion for the mem- 
brane, we still need to express the impermeability condi- 
tion 



dh q 



0,t) 



f^(M), (A3) 



where the velocity u z (k, t) is given in Eq. (|A1(I . Perform- 
ing the latter integral, we obtain the Langevin equation 
describing the relaxation dynamics of a membrane 



dh i h , f 
— - uj q h q + f q 



(A4) 



with uj q = E q /(Anq). The white noise f q is re- 
lated to the stochastic variable g(k, t) through f q (t) = 
(27r) _1 Jdk z g z (c[, k z , t). Integration of Eq. (|A2|) shows 
that its distribution is still Gaussian with zero mean 
value, (fq(t)} = 0, and correlations given by 



k B T 
2nq 
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Note that these correlations have been obtained after di- 
rect integration of the hydrodynamic equations, in the 
same way as the Langevin equation of a Brownian par- 
ticle can be found from fluctuating hydrodynamics |25| . 
In particular, the relaxation of a membrane undulation 
is characterized by (\h q (t)\ 2 ) — {\h q (0)\ 2 ) exp[— uj q t], with 
equilibrium fluctuation amplitude 
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k R T 
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which is expected from the fluctuation-dissipation theo- 
rem. 



APPENDIX B: EVALUATION OF THE GREEN'S 
FUNCTION 



In this appendix, we want to derive the correction to 
the Oseen tensor coming from hydrodynamic interactions 
with the membrane. Because the Navier-Stokes equation 
is linear, the general solution to an arbitrary force field 
f is given by Eq. In this problem, the embedding 

medium is anisotropic so that the Oseen tensor Q is not 
translationaly invariant in space. The Green function is 
obtained as the response to a point-force applied at point 
r' at instant t, f (r, t) — f (t)S(r — r'). After Fourier anal- 
ysis, the solution of the Navier-Stokes equation together 
with the incompressibility condition reads 



u(k,w) = 



nk 2 



X-kk 



{f(o 



*(k,c^)} 



( B1 ) 

with the same notation as in Appendix lAl Yet, the de- 
formation field is coupled to the fluid velocity through 
the impermeability condition Eq. (|1(J|) . This condition 
is more conviniently expressed as iujh q {uj) = u z (q, 0,to), 
where we used the set of variables {q, z, uj}. Finally, com- 
ing back to the variables {k, lu}, we find the following 
expression for the restoring force density 



(k,w) 



.E a 



dk z 
2tt 



u z (k,w) 
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Obviously, this result together with Eq. (|B1|) gives an in- 
tegral relation for the velocity field. This can be solved 
by writing explicitcly the equation satisfied by the z- 
component of the velocity. Integrating both side of the 
equation over k z , we get 



J 



dk z . 
2^ 



(k,w) 
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with oj q = E q /(4r/q) and p' = (x',y r ). Finally, bringing 
together the results Eqs. (fB If) . (|B2|I and (|B3(I shows that 
the velocity field is proportional to the force, Uj(k, u>) = 
^■5ij(k,w)/j(w). The Oseen tensor can then be writ- 



ten as Q = Qq + with the isotropic contribution 
£o(k, oj) = (X — kk)/(^/c 2 ). The explicit expression for 
the second contribution is given in Eq. (|12fl in terms of 
the variables {q, z, uo}. 
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